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Abstract:  This  report  discusses  a  Cranl'C-Nlcolson  type  of  finite 
difference  formulation  for  the  problem  of  transient  temperature 
distribution  In  a  uniform  cylindrical  pipe  with  external  Insulation 
and  containing  a  fluid  for  which  temperature,  pressure,  and  flow 
rate  are  given  as  functions  of  time.  All  material  properties  are 
assumed  to  depend  on  temperature.  Discussion  Is  also  given  of  a 
digital  computer  program  v;hich  Implements  this  formulation.  This 
implementation  incorporates  properties  of  superheated  steam,  the 
internal  fluid,  and  of  2-1/4  Cr  1  Mo  (P22),  the  pipe  material. 
Stress  calculations  are  also  treated. 
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1.  Objective 

The  purpose  of  this  monograph  Is  to  present  and  to  document  a 
procedure  for  detennlnlng  transient  temperature  distribution  and  the 
corresponding  significant  thermal  stress  components  in  metallic  pipes 
which  are  insulated  externally  and  v/hlch  contain  steam  the  temnera- 
ture,  pressure,  and  flow  history  of  which  is  specified.  Althoup^  a 
classical  solution  is  known  (13)*  for  an  Idealized  version  of  this 
problem,  it  is  not  canable  of  dealing  with  the  variable  heat  trans¬ 
fer  coefficient  which  governs  the  exchange  of  heat  between  the  nipe 
wall  and  the  flowing  steam  it  contains;  this  coefficient  is  a  comp¬ 
licated  function  of  the  properties  of  steam.  Furthermore,  the  an¬ 
alytic  solution  presumes  constant  metal  thermal  conductivity  and 
thermal  diffusivlty  while  actual  piping  materials  exhibit  very  sig¬ 
nificant  variations  of  conductivity  and  diffusivlty  over  ranges  of 
teirperature  usually  encountered  in  engineering  practice.  Accordingly 
a  numerical  arjproach,  capable  of  dealing  with  these  variations,  has 
been  chosen  for  the  procedure. 

2.  Derivation  of  the  governing  nonlinear  partial  differential  equation 

A  convenient  point  of  departure  for  the  derivation  of  the  govern¬ 
ing  nonlinear  partial  differential  equation  is  the  relation  (  5  ) 


,  9T  9  /jjy.  9T,,  ^  9  9T,.  ^  9  9T,. 


(1) 


Here  and  later  T  denotes  teirperature,  t  denotes  time,  c  denotes  spec- 


*  References  are  listed  on  pages  28  and  29. 
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ific  heat  at  constant  pressure,  and  p  denotes  mass  per  unit  volume. 

Here  C,  n,  and  c  denote  ReneralizecJ  space  coordinates  and  p,  and 

V  are  coefficients  such  that  the  distance  ds  betv/een  points  (f, ,n,r,) 

and  ( ^+dC,n+dn ,C+dc)  is  given  by 

(ds)^=  X^(d5)^  +  y^(dn)^  +  v^(drj^  (?) 

Specializing  to  the  cylindrical  coordinate  system,  v;e  take  C  =  r, 

n  =  B,  ^  =  z.  liiirthernore ,  we  assume  axial  symmetr:/  so  that  -^  =  n, 

and  we  presume  no  variation  with  respect  to  axial  position  so  that 
9T 

=  0.  V/e  thus  easily  arrive  at  the  evaluations  A  =  v  =  1,  ij  =  r. 


and  at  the  equation 
1  8T  _ 


_ _  +liT  iak^ 

a  8t  9r^  r  9r  k  9r  9r 


(3) 


vjhere  a  =  k/pc  is  the  diffusivity.  The  variation  of  k  with  respect 
to  r  is  due  solely  to  the  fact  that  conductivity  is  a  function  of 
temperature  which  itself  depends  on  radius,  'fhus  we  finally  arrive 

1  9T  _  d^T  ^  1  9T  ^  ,k',,9T  2 
a  9t  9i^2  p  ^'9r 

where  k*  =  dk/dT  for  the  pipe  material.  If  the  variation  of  k  with 
respect  to  T  were  simple,  there  would  be  some  point  in  dealing  with 
the  Kirchhoff  transformation  (6).  However,  it  is  not,  and  we  will 
deal  directly  with  the  nonlinear  partial  differential  equation 
3.  Finite  difference  formulation 

V/e  employ  a  Crank-!  Ji  cols  on  finite  difference  f'onnulation,  as 


follov/s.  Equation  ^  is  differenced  in  time  according  to 
[(HT)"^  +  (1IT)]/2  =  (T'^  -  T)/[6(rt'^  +  rt)/2] 


(5) 


where  6  represents  the  time  increment,  superscript  (  )  indicates 
evaluation  at  the  advanced  time  t  +  6 ,  and  H  is  the  space  operation 
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(6) 


The  latter  is  next  represented  approximately  by  a  finite  difference 
operation.  Vfe  use  the  notations 


9^T 

A  =  radius  increment,  A  =  — 5- 

Sr'^ 

(Note  the  new  use  of  symbol  p.)  Then^  CGi^erally, 

+  BA  +  AAV2;  -BA  +  AAV2  (7) 

where  subscript  (^)  indicates  evaluation  at  r  =  r^  =  n^  ,  subscript 
indicates  evaluation  at  r  =  r^_^  =  (n^-1)  ,  etc.  (In  general  n^  will  not 
be  an  integer.)  We  easily  find 

A  =  (T^^^  +  -2T^)/a2;  B  =  (T^^^  -  T^_^)/2A  (B) 

so  that 

Setting  6^  =  4A^/6(at+a^) ,  we  arrange  the  general  equation  in  the  form 

[-1  +  —  +  pt(T'I‘.T-Tt  J/iljTt  ,  +  (6.+2)T'I' 

2nj^  1  1+1  1-1  1-1  1  1 


o  _  9T 
‘  ^  "  9? 


,  p=  k'A 


(9) 


i+1 


=  (8,-2)T,  + 

On  the  left,  the  bracketed  coefficients  of  the  advanced  tenperatures 

T^_]_  Tt_j_^  themselves  contain  these  (unknown)  advanced  temperatures. 

Accordingly,  it  is  contenplated  that  iteration  must  be  used  in  which  earlier 

“f* 

evaluations  of  and  are  used  in  evaluating  the  bracketed  coeffi¬ 

cients  . 

Equation  10  is  valid  for  Interior  nodal  points,  i  =  2,3,...,N,  where 
N  is  the  nuirber  of  subdivisions  of  the  pipe  wall.  At  the  inside  surface 
where  the  radius  is  r^  =  n^A,  there  is  a  convective  boundary  condition 


-I4- 


(11) 


Where  (p  Is  the  fluid  tenperature  and  h  is  the  surface  heat  transfer 
coefficient.  We  vfrite  this  aa 

=  B  =  (o/A)(T  -(^)  (12) 

where  a  =  hAA. ,  and  we  also  use  the  first  of  equations  7,  takinp; 
i  =  1.  Ihus  we  get 

=  {2(T2-T,)  -  a(T,-(t))[2  -  ^  -  Dja(Tj-(l))  D/A^  (I3) 

Thus,  the  appropriate  equation  for  the  inner  surface  is 
{3i  +  2  +  a‘^[2  -  ~  -  o‘^a'^(T'[’-2<>‘^)]}Tt  -  2T2 

=  (3,-2)Ti  +  2T2  -  o(Ti-(}))C2  -  ~  -  D,o(T  -6)] 

Ui  1 

+  (2  ~  ^  +  p'^a''"(()^)  (1^1) 

where  3i  =  4A^/6(a'|^+aj ) .  Again,  note  the  nonlinearity  represented 
by  the  appearance  of  the  advanced  temperatures  in  the  bracketed  coef¬ 
ficients  and  on  the  rif^t  side.  The  advanced  tenperature  also  appears 
in  Bj  and  in  since  and  h  depend  on  the  advanced  inside  surface 

metal  tenperature. 

Finally,  the  last  equation  pertains  to  the  exterior  surface  node 

9T 

(numbered  M  =  N+1)  where  there  is  perfect  insulation  so  that  “  0* 

Employing  also  the  second  of  equations  7  with  i  =  M,  we  get 

UT  =  2(T^-Tj^)/a2 

and  the  last  difference  equation  becomes 

+  2T„ 

A.  flethod  of  solution 

Because  of  the  nonlinearities,  an  iterative  solution  is  used. 
However,  since  the  coefficient  matrix  is  tridiagonal,  it  is  convenient 
and  expeditious  to  use  a  pivotal  method  several  times  rather  than  to 


(15) 

(16) 


ar 
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use  a  purely  iterative  method  such  as  that  of  Gauss-.Seidel.  Vfe  mav 
think  of  the  equations  in  the  form  ( cf .  Reference  l^t  ) 


^  ^  .  dj  (1=1,2,.. 

(17) 

with  a^  =  Cj^  =  0. 

Then  we  form 

o 

11 

O 

c|  =  Cj_/(b^-a^c|_^)  (i=2,...,M) 

(18a) 

O 

II 

df  ^^l"^i^i-l^/^^i"^\^i-l^  (1=1,2,. 

. .  ,M) 

(I8b) 

and  obtain  the  advanced  tenperatures  as 

II 

tN  d|  -  c|T^^^  ( 1=N  ,N-1  ,...,2.1) 

(19) 

Because  of  the  appearance  of  the  advanced  temperatures  in  the  coef¬ 
ficients,  several  Iterations  are  required. 

5.  DlRital  corrputer  program. 


The  digital  conputer  program  discussed  in  appendix  C  hereto 
implements  the  theory  developed  above.  The  main  program  is  called 
PIPETEM.  The  properties  of  the  fluid  are  determined  by  subroutines 
CSLBP,  PRAND,  and  ABSW.  These  are  called  by  subroutine  AITCH  which 
determines  the  surface  heat  transfer  coefficient  h.  The  only  fluid 
whose  properties  are  incorporated  in  the  present  program  is  superheated 
steam.  The  properties  of  the  metallic  piping  material  are  determined 
by  subroutines  METAL,  CDEIXP,  and  BYUNG.  All  these  subroutines  will 
be  discussed  further  in  appendix  A  hereto. 

6.  Properties  of  the  solution. 

Having  determined  the  advanced  temperatures  T^,  (i=l,...,M), 
certain  significant  properties  of  this  solution  are  determined  in  a 
subroutine  called  AVEIU.  These  properties  are  TAV,  DTI,  ITTP,  SL, 
and  SR  which  will  be  defined  and  described  at  this  point. 
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TAV,  OTl,  and  OT2  are  rigorous  counteroarts  to  "average  tenper- 
ature,"  "delta-tee-one,"  and  "delta-tee-two"  given  in  various  pining 
and  vessel  codes  for  the  estimation  of  thermal  stress.  As  will  he 
pointed  out  in  appendix  B,  for  uniform  cylinders  with  axial  symmetry 
of  loading  and  temperature,  the  proper  stresses  to  be  used  in  code 
evaluations  may  be  derennined  directly  without  the  use  of  these  quan¬ 
tities;  however,  their  use  is  widespread  and  they  are  calculated  for 
conparison  purposes.  The  present  writer  is  responsible  for  the  defi¬ 
nitions  given  in  the  codes.  For  simplicity,  these  definitions  were 
based  on  flat  plate  geometry  and  this  is  quite  reasonable  for  pipes 
except  those  having  unusuall'/  thick  walls.  That  is,  in  these  defin¬ 
itions,  pipe  diameter  was  assumed  to  be  quite  large  comparted  to  nine 
wall  thickness.  In  subroutine  AVERT  this  assurption  is  not  made.  The 
appropriate  definitions ,  based  on  truly  cylindrical  geometry ,  are  given 
below,  followed  by  demonstrations  that  they  reduce  to  the  code  defin¬ 
itions  as  pipe  radius  goes  to  infinity  while  malntailning  constant  wall 
thickness . 

TAV,  also  designated  as  is  the  average  metal  temperature  de¬ 


fined  as 


as 

^b 

’b  * 

T  =  /TdA//dA  =  2tt 

Trdr/27T 

rdr  =  (1/hc) 

a 

a 

fb 


Trdr 


(20) 


where  a  =  inside  radius,  b  =  outside  radius,  h  =  b-a  =  wall  thickness 
(do  not  confuse  with  earlier  use  of  h  as  heat  transfer  coefficient). 


and  c  =  (b+a)/2  =  mean  radius.  If  we  define  y  =  r-c,  then 
|■+h/2  C+h/2 

T  =  (1/h)  (l+y/c)Tdy  ->  (1/h) 

‘’-h/2  ^-h/2 

and  this  is  the  definition  given  in  the  codes. 


My  as  c  -*■  ” 


(21) 
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V/e  next  define  a  quantity  which  can  be  tenned  the  thermal 


+h/2 


moment ,  viz . : 

M  =  fr(T-T)dA  =  Pirf  i^T-T)dr  =  (l+y/c)==(T-T)dv 

^  J  a  J -h/2 

Using  the  definition  of  T,  this  can  be  simnlifled  to  the  form 


M  =  2ttc 


r+h/2 

( 

-h/2 


(l+y/c)'IVdy  -  7TThV6 


(22) 


(23) 


Tne  next  step  in  this  development  is  to  consider  the  average  and 
the  thermal  moment  of  a  linearly  varying  temperature  distribution 
which  has  variation  V  from  r  =  a  to  r  =  b.  Such  a  temperature  is 
given  by  T  =  Vr/h  for  which  we  find 


T  =  (V/h^c) 


fb 

r^dr  = 


(V/3h^c)(b^-a^) 


(2/1) 


fb 

M  =  2ti  r^  ( Vr/n-T)  dr  =  ( TtVch  V6 )  ( l-h^  /12c  ^ )  ( 2S ) 

■'  a 

It  requires  considerable  manipulation  to  obtain  the  last  form. 

By  equating  these  two  values  of  M,  one  evaluates  the  (equiva¬ 
lent)  linear  variation 


V  =  [12 


;+h/2 

( 

-h/2 


( l+y/c)'IVdy  -  ThVc]/h^(l-h^/12c^) 


f+h/2 

(12Ai^)  Tyd:/  (26) 


^-h/2 

as  c  The  limiting  values,  as  c  are,  in  each  case,  the 

definitions  given  by  the  various  piping  and  vessel  codes.  Ttie  present 
subroutine,  AVERT,  enploys  the  more  general  definitions  based  on 
cylindrical  goonetry  with  finite  radii.  The  variation  V  is  also 
denotes  as  ATj  in  the  codes  and  as  CTl  in  the  output  of  the  present 
program. 
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In  the  codes  the  quantity  AT^  is  defined  as  the  nonlinear 
residual  terperature  at  the  inner  or  the  outer  radius  (v/hichever 
ftives  the  greater  result)  after  the  constant  part,  T,  and  the  linearly 
varying  part,  V,  are  subtracted.  In  the  outnut  of  the  present  program 
this  quantity  is  denoted  as  1712  and  defined  as 

012  =  rtoc{  1  Tj^-T-Ori/2 1  , 1  T^-T+Ori/2 1  }  ( 27 ) 

The  quantities  ATj  and  AT^  are  Intended,  in  the  codes,  to  permit 
estimating  significant  thermal  stresses  in  general  pressure  vessel 
components.  However,  in  the  particular  case  of  uniform  hollo.'/  cylin¬ 
ders  ,  the  stresses  resulting  from  radial  temperature  gradient  are 
definitely  classified  as  "local  thermal  stresses"  and  there  is  no 
actual  need  to  separate  the  actual  temperature  distribution  as  has 
been  indicated  above.  Instead,  the  thermal  stresses  can  be  calculated 
directly,  and  this  is  done  in  subroutine  A''/ERT.  The  quantities  SL 
and  SR  are  the  values  of  circumferential  (or  axial  —  they  are  enual) 
stress  due  to  the  temperature  distribution  at  inside  and  outside  rmtal 
surfaces ,  respectively ,  under  the  assumption  of  zero  net  axial  force 
and  bending  and  twisting  moment.  The  latter  contributions  may  be  de¬ 
termined  separately  by  use  of  a  conventional  pining  flexibility  analysis. 

Because  of  the  variability  (with  temoerature)  of  the  Young’s 
modulas  and  of  the  coefficient  of  thermal  expansion,  the  calculations 
are  not  trivial.  Discussion  of  this  matter  is  given  in  appendix  B 
hereto. 
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APPEI'JDIX  A  Discussion  of  the  subroutines 


Subroutines  CSUBP,  PFJAND,  and  ADSIIU  receive  as  inputs  the  tem¬ 
perature  (°F)  and  the  pressure  (psia)  of  the  steam  flcAfinp,  In  the  pine 
and  produce  as  outputs  the  values  of  Cp  (specific  heat  at  constant  pres¬ 
sure),  Up  (Prandtl's  number),  and  p  (absolute  viscosity)  or  the  stefn. 

The  calculations  are  by  polynomials  in  T  and  P  the  coefficients  of  which 
were  obtained  by  least  squares  analysis  of  data  points  read  from  the 
graphs  on  pages  293,  29^,  and  297  of  the  ASMB  Steam  Tables  (  ).  It  is 

inpllcltely  assumed  that  the  steam  is  saturated  or  superheated.  If  the 
steam  is  wet,  the  proper  values  are  probably  not  obtained. 

Sixty  number  pairs  were  used  in  dtermining  the  coefficients  in 
CSIBP,  fifty-four  pairs  for  PRAUD,  and  twenty-four  pairs  for  ABS'iU.  The 
only  type  of  accuracy  evaluation  perfonned  v;as  to  observe  that  each  or 
the  "subroutines  was  able  to  generate  the  input  numbei’  pairs  within  one 
or  two  percent  (generally  much  closer  than  this).  Such  accuracy  is 
greater  than  that  implicit  in  convective  heat  tranfer  theory  and  certainly 
at  least  as  great  as  the  accuracy  with  which  the  input  values  were  read 
from  the  graphs. 

Subroutine  AITCH  determines  the  film  heat  transfer  coefficient  h 
by  use  of  the  Colburn  formula  (cf.  Gledt  (  8  )) 

h  =  Ug^GCp  (A-1) 

where  G  is  mass  per  unit  time  per  unit  area,  i.e.. 


G 


=  m/A 


flow 


(A-2) 


and  Uom  is  the  Stanton  number, 

oi 
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(A-3) 


N  =  0  02^  rj~2/3 
ST  'TR 

In  thiy  formula  is  the  Prandtl  nurtjer  and  Is  the  Reynolds  number 

where  d  is  the  pipe  inside  diameter.  The  Colburn  fomula  is  one  which 

is  recommended  for  cases  of  well  develoned  turbulent  flow  at  fairly 

hipji  Reynolds  nuirhers.  It  seems  to  be  an  apnropriate  choice  excent 

perhaps  for  cases  in  which  the  operatinp;  transient  involves  a  ver'/ 

{■pi^at  decrease  in  mass  flow  rate.  In  such  cases  Q  indeed  aporoaches  zero 

but  not  necessarily  in  the  way  reoresented  by  the  Colburn  formula. 

Incidentally,  subroutine  AITQi  avoids  a  comoutational  difficulty  as 

G  goes  to  zero  by  corrhining  the  equations  in  the  form 

h  =  0.023  Cp  G°*^  (A-5) 

The  specific  heat  c  is  evaluated  at  the  fluid  bulk  temoerature 

P 

but  rjpp  and  p  are  evaluated  at  the  "film"  temperature  which  is  the 
arithmetic  average  of  fluid  bulk  temperature  and  inside  wall  metal 
tenperature.  Since  the  latter  is  an  unlmo’wn  in  the  calculations,  an 
iterative  evaluation  must  be  used. 

Subroutines  rCTAL,  COEXP,  and  EYUIG  obtain  pertinent  properties 
of  the  metal  pipe  material  corresponding  to  input  values  of  the  metal 
temperature.  These  subroutines  presently  provide  data  for  only  one 
particular  material,  the  standard  low  chrome-mo ly  hiph  temperature 
alloy  designated  by  the  number  22  —  pipe  is  deslmated  as  P22.  Sub¬ 
routine  METAL  provides  values  of  thermal  conductivity  (k),  the  rate 
of  variation  (k*  =  dk/dT),  and  the  thermal  diffuslvlty  using  poly¬ 
nomial  representations  of  material  test  data. 
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Althou^  Table  1-4.0  of  the  Nuclear  Components  Code  (2) 
sup;;ests  using  the  same  values  for  thermal  conductivity  and  thermal 
diffusivity  for  P22  as  for  low  carbon  steel,  data  given  in  the  TT’RC 
Data  series,  References  (15)  and  (16) »  indicate  a  sufficiently  great 
dependence  on  alloying  constituents  that  we  have  preferred  to  use  the 
latter.  Thus  we  get  substantiallv  lower  values  of  conductivity  and 
diffusivity  than  the  ASflE  tabulation  gives,  '[he  actual  data  employed 
in  constructing  subroutine  F'ETAL  are  the  data  points  from  curve  19 
page  1158  of  Ref.  (15)  and  curve  2  page  343  of  Ref.  (I6).  It  should 
be  noted,  however,  that  it  v/ould  be  a  very  sinple  matter  to  modify 
the  subroutine  ICTAL  to  use  the  ASFC  data  rather  than  the  TPRC  data. 

Subroutines  CCEXP  and  EYUNG  are  used  (by  subroutine  AVB^RG)  to 
provide  the  coefficient  of  thermal  expansion  and  the  Young's  modulus 
corresponding  to  the  given  input  temperature.  The  data  is  from 
Appendices  D  and  C  of  the  Power  Piping  Code,  Ref.  (  3),  and  linear 
interpolation  is  used.  Poisson's  ratio  is  assumed  to  have  the 
constant  value  0.3. 

All  subroutines  as  well  as  the  main  program  use  English  custom¬ 
ary  units,  specifically  Inches  (rather  than  feet)  and  seconds  (rather 
than  minutes  or  hours). 
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APPEMDIX  B  Stress  Calculations 


Assuming  what  is  called  a  "quasi-static"  situation,  i.e., 
neglecting  forces  associated  with  accelerating  the  mass  of  the  pipe 
as  its  thermal  expansion  causes  displacerrent  of  its  material  particles, 
and  assuming  also  that  the  difference  betv/een  adiabatic  and  isothermal 
material  constants  is  of  neglegible  significance,  vje  nevertheless  are 
faced  with  a  nontrivial  problem  in  determining  the  state  of  stress 
which  results  from  the  terrperature  changes  discussed  in  the  body  of 
this  report. 

Two  material  constants,  E  (Young's  modulus)  and  ot(the  coefficient 
of  linear  expansion  —  do  not  confuse  with  thermal  dlffusivity)  are  of 
significance  in  connection  with  stress  evaluation.  Both  vary  with 
terrperature.  Further,  as  will  be  discussed  below,  ot  varies  with  stress 
state.  A  number  of  writers  have  dealt  with  stresses  in  tubes  due  to 
chanp^s  in  terrperature,  among  them  Hilton  (  9)  and  Chang  and  Chu  (7)* 


The  analysis,  in  effect,  reduces  to  the  following. 

=  [E/r2(l-v)][(r2-a")  J(b)/(b2-a2)  -  J(r)]  (B-1) 

=  [E/r2(l-v)][(r2+a2)  J(b)/(b2-a")  +  J(r)  -  r^oTCr)]  (B-2) 

U 

0  =  Ee  +  [E/(l-v)][2vJ(b)/(b^-a^)  -  ofr(r)]  (B-3) 

z  z 

where 

J(r)  =  [  CaT(C)dC  (B-'O 

a 

Presuming  that  the  axial  force  P  is  zero,  we  get 

G  =  2  J(b)/(b2-a^)  (B-5) 

z 


and  can  thus  evaluate  o  .  For  all  temperature  distributions  which 

z 
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result  from  external  heating  (i.e. ,  from  heat  applied  to  the  Inner  or 
outer  surfaces  but  not  generated  within  the  wall  Itself)  the  extremal 
stresses  v;lll  be  at  r  =  a  or  at  r  =  b.  (With  external  Insulation,  It 
Is  plausible  that  the  extreme  must  be  at  r  =  a.)  Vfe  thus  have  =  n 
at  r  =  a  and  at  r  =  b,  and  we  also  have 

=  [E/(l-v)][2J(b)/(b^-a^)  -  ntT(a)]  at  r  =  a  (B-6) 
Og  =  =  [E/(l-v)][2J(b)/(b^-a^ )  -  rtT(b)]  at  r  =  b  (B-7) 

In  computing  J(r)  and  J(b),  the  a  under  the  sign  of  Integration 
Is  the  "Instantaneous"  value  appropriate  toT(r).  It  Is  (nresumably) 
the  "zero-stress"  value  of  a  which  Is  tabulated.  Obviously,  hovjever, 
during  the  temperature  changes  the  various  elements  of  the  pine  wall 
are  not  stress-free  and  the  question  arises  oO  hoi^  It  Is  that  we 
should  use  the  zero-stress  value  of  a.  The  explanation  is  surety  not 
well  known  (this  question  is  not  even  considered  in  any  of  the  standard 
references  on  thermal  stress  analysis);  the  explanation  v/hich  follows 
was  given  by  the  writer  in  a  report,  dated  Noveirber  i960,  entitled 
"Seflections  on  Piping  Flexibility  Analysis,"  v/hich  was  intended  to  ex-' 
plain  the  basis  of  procedures  incorporated  in  the  "MDC  Docianent"  of 
the  rAechanlcal  Design  Committee  of  the  ANSI  (then  ASA)  B31  Piping  Code 
Conmittee. 

This  explanation  considers  a  one  dimensional  constrained  thermal 
expansion  problem  of  a  carbon  steel  bar  of  length  L  which  at  70  ^ 
fits  precisely  between  two  rigid  and  immovable  anchors;  that  is,  at  70  F 
there  is  no  stress  and  no  gap.  The  temperature  is  increased  to  300  P. 
One  is  to  determine  the  resulting  compression  stress,  presuming  that 
there  is  no  tendency  to  buckle  laterally.  As  fundamental  we  take  the 
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relation 


de  =  da  a:  +  adP 


(15-8) 


and,  since,  in  this  problem  de  =  0,  v;e  have 

da  =  -EadT  (15-9) 


a 


P 


f300 

E(T)a(T)dT 

'70 


vjhere  subscript  (p)  indicates  the  "final"  condition. 


(F5-in) 


Appendices  B  and  C  of  AfISI  B31. 1-1973,  Ref.  (  ),  r,ive  the 

follavlnR  data  for  carbon  steel  pipe: 


T 

E‘10"^ 

a*  10^ 

Ea 

70 

27. 

6.07 

169.35 

200 

27.7 

6.38 

176.73 

300 

27. '1 

6.6o 

180.8^ 

Using  trapezoidal  Integration,  we  evaluate  equation  (B-10)  to  get 


bp  =  -(169. 35+176. 73)  (130/2)+(  176.73+180.8^1)  (100/2) 


=  -11037^  psi  (B-11) 

(This  result  is  incorrect  as  we  will  show.) 

An  alternate  calculation  obtains  the  answer  by  considering  one  anchor 

to  be  removed  and  letting  the  bar  expand  freely  (fST^-stress  condition) 

from  70  P  to  300  P  and  then  supplying  sufficient  compressive  load  to 

restore  the  bar  to  its  original  length.  This  calculation  gives 
^  (300  . 
e  =  (T)dT  =  [(6.07+6.3B)(65)+(6.38+6.60)(50)]-10"^ 

■'70 

=  1. '4583*  10"^  (B-12) 

Op  =  -(27.^*10^)(1.'45B3*10“^)  =  -39956  psi  (B-13) 

which  is  a  slightly,  but  definitely,  smaller  value. 

Now  it  would  Indeed  be  disturbing'  if,  without  any  inelastic 

behavior,  the  final  stress  state  were  to  denend  on  the  sequence  of 
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load-temperature  application.  This  would  indeed  be  a  new  type  of 
thennal  ratcheting.  However,  we  will  new  shm  that  the  second  cal¬ 
culation  gives  the  correct  result  for  any  load-tenperature  sequence, 
based  upon  the  postulate  that  the  final  stress  state  is  independednt 
of  sequence  and  depends  only  on  material  properties,  the  final  tem¬ 
perature  state,  and  the  final  strain.  In  other  words,  wc  presume  that 
equation  (B-  B)  represents  an  exact  differential  so  that 

Sp  =  fCOpjTp) 


If  this  is  the  case,  we  have 


3  /  3e \  _  3  /  3e\ 
It  da’  ^  ■^3T^ 


and,  since 


we  have 


3e 

3T  “ 


^  =  i 

3a  E 


(B-li4) 


(B-15a,b) 


3a  _  d  /1-, 
3a  dT'^E'^ 


(B-16) 

(v;e  have  written  d  rather  than  3  on  the  ripht  since,  by  its  definition, 
E  does  not  depend  on  a . ) 

Ihus,  if  E  varies  with  T,  then  a  must  vary  with  a;  i.e.,  a  = 
a(T.a).  The  tabulated  values  we  use  are  values  of  a(T,0)  and  these 
were  appropriate  in  the  second  numerical  calculation  but  not  in  the 
first.  Thus,  for  our  original  problem  we  should  integrate 

da  =  -E(T)o(T,a)dT  (B-17) 

taking  into  account  the  dependence  of  o  upon  a.  V/e  have 

a(T,a)  =  a(T,0)  +  =  a(T,0)  +  a^i)  (B-lB) 

so  that 


da  =  -E(T)[a(T,0)  +  (a/E^)(dE/dT)]dT 
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=  -r':(T)a(T,0)dT  -  oclE/E 


(B-19) 


Thus , 

-a(T,0)dT  =  (Eda  -  adE)/^’^  =  d(o/E)  (n-20) 

Performing  this  integration,  v;e  get 

f300 

(o„/E„)  =  -  n(T,0)dT  = 

■'70 

Op  =  (-1JI583‘10"^)(27JM0^)  =  -39956  psi 
as  in  the  second  calculation. 

In  the  more  conplicated,  multi-axial  stress  situation  with  which 
this  monograph  is  primarily  concerned,  the  argument  takes  the  followinf’; 
form.  We  Imagine  any  element  (dr,rd  ,dz)  of  the  pipe  to  be  separated 
from  the  adjacent  elements  and  to  exoerience  stress-free  thermal  ex¬ 
pansion  resulting  in  extensional  strains  of  the  form 

fTa 

e  =  a(T,0)dT  (D21) 

h: 

Then  surface  forces  are  anplled  to  the  element  to  cause  it  to  assume 
a  form  sudi  that  all  elements  can  be  fitted  to  gather  in  the  final 
stressed  state.  Since  equations  (D-l)  throuph  (B-7)  were  derived  by 
satisfying  the  equilibrium  and  constitutive  laws,  they  indeed  describe 
the  final  state. 

One  final  remark  needs  to  be  made  concerning  stress  calculations 
for  code  purposes.  Althou^  the  so  called  "simplified"  methods  of  an¬ 
alysis  provide  for  stress  comoonents  corresponding  to  the  temperature- 
like  quantities  ATj  and  ATj ,  it  is  very  important  to  note  that  in  the 
specific  case  of  a  uniform  pipe  or  cylindrical  vessel  the  stress  state 
resulting  from  a  radial  tenperature  gradient  does  not  have  to  be  senar- 
ated  into  "secondary"  and  "local"  parts  but.  Instead,  the  actual  stress 
may  all  be  considered  as  "local."  The  v/ord  "local"  is  quite  imoortant 
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in  considering  allov/able  states  of  stress,  A  definitive  staterrent  to 
this  effect  is  to  be  found  on  the  fifth  page  of  the  AS'E  "Criteria" 
document ,  Ref.  (  *  viz : 

"A  special  exception  to  these  general  rules  is  the  case  of  the 
stress  due  to  a  radial  tenperature  gradient  in  a  cylindrical 
shell.  It  is  specifically  stated  in  II-^(12(m)(2)(6)  of  Section 
III,  and  in  4-112(1) (2) (6)  of  Appendix  4  of  Division  2,  that 
this  stress  may  be  considered  a  local  thermal  stress.  In 
reality  the  linear  portion  of  this  gradient  can  cause  defor¬ 
mation,  but  it  was  the  opinion  of  the  Special  Committee  that 
this  exception  could  be  safely  made." 

'fhus  the  stresses  calculated  according  to  equations  (B-6)  and  (I.i-7) 
of  this  appendix  themselves  directly  give  the  significant  thermal 
stresses  due  to  the  radial  gradient,  and  these  are  indeed  local 
thermal  stresses. 
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APPEfJDIX  C  Instmctions ,  Data  Preparation .  Tynical  Problem 


Instructions  for  the  use  of  the  programs  and  preparation  of  in¬ 
put  data  are  given,  in  the  form  of  comments,  in  the  initial  portion 
of  the  program  Itself.  This  portion  of  the  program  reads  as  follOTs: 
Temperature  distribution,  etc.,  in  uniform  pipe  with  exterior 
insulation,  convective!:/  heated  by  fluid  contents.  The  basic 
procedure  is  a  flnte  difference  treatrrent  of  the  one  (space) 
dimensional  partial  differential  equation  in  cylindrical  co¬ 
ordinates.  r-letal  properties  are  assured  to  vary  witli  temper¬ 
ature.  Surface  heat  transfer  coefficient  is  determined  by 
Colburn’s  relationship.  The  fluid  content  in  this  version  is 
assumed  to  be  superiieated  steam  and  the  pipe  material  is  P22. 

Input  is  by  neans  of  two  kinds  of  cards.  The  first  card  reads 
OD,  liTT,  TINIT,  N,  NIT,  and  HPROB  according  to  FORriAT(3FiO.0 ,3IR) . 
OD  is  the  outside  diameter  in  inches.  VJT  is  the  wall  thickness 
in  inches.  TINIT  is  the  uniform  initial  metal  temperature  in 
degrees  Fahrenheit.  (For  nonunifom  initial  tenperature,  see 
below.)  N  is  the  (even)  number  of  equal  subregions  into  which 
the  wall  thickness  is  divided.  (If  input  N  is  odd,  it  is  round¬ 
ed  up  to  the  next  even  number.)  NIT  is  the  nunber  of  iterations 
errployed  each  time  step  in  order  to  solve  the  nonlinear  matrix 
equation.  fJPFiOB  is  an  integer  identification  number. 

If  TINIT  <  -500.,  the  program  looks  for  nonuniform  initial  metal 
terrperature ,  reading  enouj^  cards  according  to  'F0Rr'lAT(7’^10.0) 
to  provide  the  (N+1)  required  initial  temperatures. 
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Subsequent  cards  read  I,  MIEES,  TIML,  TIM2,PHI1,  PHI2,  PRESl, 
PRES2,  FLOWl,  and  FI£)W2  according  to  F0RriAT(2I5,8P8.0) .  I  is 
an  identif>/ing  serial  niurber  which  is  useful  for  orderinr  the 
data  but  which  is  not  used  in  the  nrogram.  MTEES  is  the  nurrber 
of  time  subintervals  into  which  the  large  time  interval  from 
time  =  TIMl  until  time  =  Tir^  is  to  bo  divided.  PHIl  and  PHI2 
are  values  of  the  fluid  temperature  at  times  Tllil  and  Tll'i? 
respectively.  This  tenperature  is  assumed  to  vary  linearly 
with  time.  Similarly  for  PRES  (fluid  pressure  in  pslg)  and 
PIDW  (fluid  flow  rate  in  pounds  per  second). 

The  data  deck  for  a  problem  terminates  with  a  card  containing 
10000  in  the  first  five  positions.  A  nuirber  of  problems,  each 
having  a  data  deck  as  described  above,  may  be  stacked.  The 
sequence  of  problems  is  normally  termnated  by  a  card  contain¬ 
ing  -1.  in  the  first  three  positions. 

Sane  thought  was  given  to  providing  an  option  which  would  permit 
several  time  intervals  to  be  calculated  per  output  line  of  print. 

It  was  decided  that  not  only  would  this  be  a  needless  complication 
but  also  that  frequent  output  is  appropriate  and  useful  during  those 
periods  when,  for  one  reason  or  another,  small  time  Increments  are 
enployed. 

It  is  easily  possible  to  make  minor  internal  changes  so  that  the 
nodal  tenpe natures ,  rather  than  the  results  calculated  in  AVERJ, 
are  printed.  This  is  useful  when  comparing  results  with  those  ob¬ 
tained  by  other  methods  of  calculation. 
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The  pror^ram  liatlng  is  not  given  here  for  reasons  of  space  econorr/. 
Tnose  who  are  seriously  interested  in  installing  the  nrogram  should 
communicate  directly  with  the  writer  to  make  appropriate  arrangements. 
The  program  is  written  in  FORTRAI'J  using  none  other  than  quite  primi¬ 
tive  capabilities.  It  has  been  run  on  an  IBn  360/67  requiring  only 
56K  core  for  execution  (both  G  and  H  compilations).  For  a  stack  o^ 
thirty  two  problems,  involving  a  total  of  b2^IO  time  steps,  with  h  = 

10  (10  layers)  and  NIT  =  5  (five  iterations),  time  requirements  \\rere 


as  follows: 

G  Compiler 

?I  Compiler 

Compile 

1^1. 38s 

23.20s 

Link-edit 

1.87s 

l,7^s 

Execute 

I4ml9.48s 

llml7.03s 

Total 

l'^m35.733 

llm^t  1.97s 
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Figure  C-1  Instnictions  for  the  use  of  PTPIiTIir'I 
'Phis  material  is  a  nart  of  the  PIPETEf'’  nromram  itself 
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APPEMDIX  D  Program  verification 


By  artificially  causing  the  program  PIPETKM  to  operate  with  constant 
values  of  diffuslvity,  conductivity,  and  surface  heat  transfer  coefficient, 
it  may  be  used  to  solve  problems  for  which  "classical"  solutions  are  known. 
Tl'ie  program  was  initially  verified  by  taking  O.D.  =  10000.,  V/.T.  =  1.,  and 
comparing  results  with  data  points  taken  from  graphs  given  by  Mcfleill  and 
Brock  (12).  vyhen  a  number  of  such  conparlsons  checked  within  the  accuracy 
with  which  the  graphs  could  be  read,  a  more  exacting  test  was  made.  This 
is  described  in  what  follows. 

Ozi^ik  (13)  gives  solutions  involving  true  cylindrical  geometry,  ^or  a 
unlfonn  cylinder  with  inside  radius  a  and  outside  radius  b,  insulated  on 
its  outer  surface,  and  convectlvely  exposed,  with  constant  surface  heat 
transfer  coefficient  h,  to  an  internal  fluid  having  "rarp"  temperature 
<j)  =  At  (where  A  is  a  constant  and  t  denotes  time) ,  and  having  initial 
tenperature  distribution  T(r,0)  =  0  at  tine  t  =  0,  the  solution  may  be 
written  as 

■Kr.t)  -  I  ^t)Vp)/A/„(l) 

m=l 

where 

T  (t)  =  2ABnt[exp(-Y_t)  -  1  +  Y  t]/a^ 

i  n  M  M 

'^ra  =  1 

B  =  ahA  (Biot's  number),  p  =  r/a,  o  =  b/a 

and  the  numbers  y  are  the  zeros  of  the  function 
m 

rBJo(u)  +  uJi(u)]Yj(pu)  -  [BYg(y)  +  yY j (y) ]J j ( oy) 

The  symbols  and  Y^  denote  Bessel  functions. 
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It  proved  to  be  quite  difficult  to  devise  a  pro^ran  to  nerform 
these  calculations.  An  apparently  different  solution,  a;iven  by  J^uikov 
(lO),  was  dealt  with  first,  but  debur/;in;^  v/as  not  conpleted  successfullv. 
A  later  work  by  Luikov  (11),  in  Russian,  nave  the  fonn  v;e  actually  used, 
citing  Ozipik  (13)  as  source.  Both  single  and  double  precision  versions 
of  an  inplementing  program  were  constructed  and  debupr,ed.  The  double 
precision  version,  much  slower  than  its  countemart,  partly  because  of 
the  use  of  a  "home-made"  DP  I3essel  function  subroutine,  was  used  to 
verify  that  the  single  precision  version  provided  adequate  accuracy. 

In  all  such  checks,  results  from  the  two  programs  agreed  to  five  sig¬ 
nificant  digits. 

The  test  problem  employed  to  check  the  integrity  and  accuracy 
of  PIPETEM  was  as  follows:  a  =  Inches,  b  =  5  inches,  A  =  1. ,  B  =  1. , 
n=  1.,  k=  1.,  k’  =0.,  h=  0.25.  Analytical  calculations  were  made 
for  t  =  0,  1.,  2.,  3.>  and  5.  PIPETEM  calculations  were  made  f’or 
t  =  o.  to  t  =  5.  at  intervals  of  0.1.  In  PIPEIEM  we  took  N  =  10  and 
HIT  =  5.  (V/e  wasted  a  few  seconds  by  not  usinf  HIT  =  1.  This  problem 
is  completely  linear  and  no  iteration  is  necessary.)  In  the  analytic 
solution  initially  we  took  account  of  ten  tems  in  the  suirmation. 

The  results  were  so  close  that  differences  can  not  be  displayed 
graphically.  However,  the  differences  were  sufficiently  great  as  to 
demand  an  explanation.  Truncation  error  in  the  analytic  solution  was 
suspected.  Accordingly,  the  analtic  solution  was  repeated  three  more 
times  using  five,  tv\fenty,  and  forty  terms.  (The  forty  term  solution 
was  very  extravagant  of  coirputer  time.)  Plotting  typical  analytic 


results  versus  reciprocal  number  of  teims  clearly  indicated  that  Indeed 
the  analytical  solution  did  suffer  from  truncation  error  and  that  extra¬ 
polation  did  indeed  lead  to  the  result  given  by  PIPETEM.  Since  the  theo¬ 
retically  "exact"  solution  has  by  far  the  greater  error ^  it  is  not  easy 
to  estimate  the  accuracy  of  the  PIPETEE  results,  but  in  this  problem  the 
error  appears  to  be  only  llttl«  n^rothan  one  part  in  one  thousand.  The 
following  tabulation  gives  calculated  results  for  T(0,5),  for  example. 

N  =  5  N  =  10  N  =  20  N  =  ^40  N  =  00  PIPETEfl 

2.0569^  2.08627  2.09991  2.106»47  2.11496  2.11220 

Table  D_i  Calculated  results  for  T(0,5) 

The  first  four  tabulated  results  ;uere  obtained  by  the  single  precision 
analytic  solution.  (The  double  precision  version  also  gave  2.0R627  for 
N  =10.)  The  extrapolated  value,  for  M  =  00  ,  was  obtained  by  assuming 
that  T  =  A  +  Bx  +  Cx^  +  Dx*,  where  x  =  5AI,  leading  to  the  evaluations: 

A  =  2.114958096,  D  =  -.073526667,  C  =  .049093333,  D  =  -.033584762.  The 
difference  between  the  extrapolated  analytic  value  and  the  value  provided 
by  PIPETEM  is  0.00275  which  is  0.13%  of  the  average  of  these  values. 

As  an  ultimate  check  on  the  validity  of  the  program,  an  artificial 
nonlinear  problem  was  made  up  and  solved  numerically.  The  conditions 
are:  inside  radius  a  =  1,  outside  radius  b  =  2,  conductivity  k  =  T  (so 
that  k'  =1),  heat  transfer  coefficient  h  =  15e^,  fluid  temperature  d 
=  6e^,  initial  temperature  distribution  T(r,0)  =  r  +  4/r.  The  diffus- 
ivity  a  is  a  function  of  both  terrperature  and  time,  viz: 

a  =  Y/[1+(4-Y)V(4+Y)^],  Y  =  (Z-/Z2-16)V4,  Z  =  Te“^ 

It  is  easy  to  verify  that  T(r,t)  =  e^(r  +  4/r)  is  a  solution  of  this 
problem.  Althouj^ti  matters  of  uniqueness  of  solution  to  nonlinear  sys¬ 
tems  are  difficult  to  deal  with,  it  happens  that  PIPETEM  indeed  yields 


the  solution  cited. 
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Attention  v/as  focussed  on  the  values  =  T(l,t*)  =  100. noo 
Tg  =  T(2,t*)  =  80.000,  where  t*  =  logg20.  Percent  errors  are 
shown  in  Table  D-2  corresponding  to  the  128  evaluations  indicated, 
values  of  r  being  1,  2;  values  of  N  being  2,  6,  10;  values  of 

NIT  being  1,  2,  5,  10;  and  values  of  DTEE  being  t*/5,  t*/20,  t*/50, 
tV200. 


DIEK  NIT  rj  =  2 


N  =  ^ 


II  =  6 


N  =  10 


tV5 


t*/20 


tV50 


tV200 


1 

2 

5 

10 

1 

2 

5 

10 

1 

2 

5 

10 

1 

2 

5 

10 


-21.51/-l'i.0iJ 
-2.726/+l.il79 
-2.592/-6.721 
-.8019/+'!.  055 
-2. 0^127+2. 726 
-1.962/-. 9932 
-1. 6^18/+.  8822 
-1.735/-. 1605 
-I.379/+2.520 
-1.8i)0/-.1536 
-1.71V+.^576 
-1. 723/+. 3710 
-I.63I/+I.020 
-I.73V+.313I 

-1. 725/+. 3715 
-1. 725/+. 3713 


-Ii4.70/-7.173 

+.3159/+3.258 

-l.7BV-6.35i 
+.5816/+3.232 
+1. 0857+^4.366 
-I.070/-I.907 
-.38W+.7397 
-.5519/-. 0690 
+.2108/+2.209 
-.7005/-. 556 8 
-.5162/+. 1779 
-.5285/+. 0778 
-.3372/+.  72  V 

-.5V6/+.023^ 

-.5287/+. 0857 

-.5288/+.04l89 


-9.532/-2.4I87 

+.0102/-. 0009 
-.9V9/-3.697 

-.02l6/+.38l4| 

+1.8440/+4|.4|25 

-1.010/-2.239 
-.1063/+. 5389 
-.2782/-. 0115 
+.6620/+2.185 
-.  W/-.5960 

-.233V+.1084I 
-.24466/-.  00  69 

-.0229/+. 6389 
-.2618/-. 0232 
-.2446V+.0365 
-.2465/+. 0351 


-3.3727+3.679 
-.72227-2. 64(6 
+I.289/+3.447 
-2.023/-6.768 
+2.5697+4. 810 
-1.104/-2.594 

+.0671/+.5541 
-.1289/-. 1902 
+.96987+2. 209 
-.3513/-. 6870 
-.0781/+. 0837 
-.0916/-. 0043 
+.1658/+. 6 182 
-.1089/-. 0567 
-.0910/+. 0131 
-.0911/+. 0114 


Table  D-2  Percent  errors  in  evaluation  of  T^^  and  T^  in  test  problem. 
Error  in  T^  (Tg)  appears  to  left  (ripjit)  of  slash  (/). 


It  may  be  noted  that  the  accuracy  is  quite  sufficient  for  engi.neering 
purposes  if  the  number  of  layers  N  =  4|^  the  number  of  iterations  NIT 
=  5>  and  the  total  time  Interval  is  divided  into  20  parts.  It  may 
also  be  remarked  that  the  total  time  to  evaluate  all  128  results  was 
only  2  minutes  20  seconds. 
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APPHJDIX  E  Heat  Transfer  at  Outer  Radius 


Pollowinfi;  the  conpletlon  of  the  precedlnp;  portions  of  this  renort, 
an  application  came  to  attention  which  called  for  the  aumentation  of 
PIPETEM  capabilities.  Ordinarily  hipb  temperature  pininp-  and  fittinns 
are  lagged  (insulated).  However,  in  some  cases  it  might  be  of  inte'^ 
est  to  follow  the  temperature  history  of  unlap-ged  pipes  v;hen  sub.iected 
to  internal  fluid  transient  temperatures.  Also,  PTPE'l’ET'I  is  obviously 
applicable  to  Insulated  flanges  if  one  discounts  axial  conduction  and 
the  additional  thermal  resistance  of  the  flange-bolt  Interfaces. 
Down-transients  (i.e.,  deceases  in  the  temperature  of  fluid  contents) 
lead  not  only  to  local  thermal  stresses  but  also  to  discontinuity 
thermal  stresses  if  there  are  significant  changes  in  ’wall  thickness, 
since  the  average  terperature  of  the  thicker  oart  lags  (in  tine)  be¬ 
hind  that  of  the  thinner  part.  This  is  particularly  tiaae  of  flange- 
pipe  assent)lies.  Leaving  the  flange  unlagged  (but  with  lapcging  on 
the  pipe  itself)  causes  the  average  temperature  of  the  flange  to  be 
lov'jer  than  that  of  the  pipe  since  the  former  can  dissipate  heat  by 
convection  and  radiation.  Then,  when  a  down- transient  occurs,  the 
resulting  thermal  dislocation  stress  will  be  less  than  v\/hat  it  would 
have  been  if  the  flange  had  been  lagged.  The  calculation  of  the 
temperatures  in  an  unlagg^ed  flange  really  calls  for  finite  element 
treatment  since  axial  conduction  and  heat  loss  ^rom  bolts,  nuts,  and 
flange  face  is  obviously  of  importance.  However,  a  useful  estimate 
of  the  effect  can  be  gained  by  considering  only  radial  flow,  with, 
perhaps,  an  artificially  enhanced  value  of  the  effective  convective 
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external  heat  transfer  cceffioient. 


Accordingly,  PIPETEH  has  been  modified  to  account  for  a  convective 
condition  at  the  outer  boundary,  viz. 

kj,^3  =  h(.!;-'j:^^)  g-l 

rather  than  B  =  0.  The  ilth  finite  difference  equation  becomes 

+  2  +  g\?.  i  -  dt;;  -  2-;; 

=  (8„-2)T„  +  2T„  -  o(T,,,-*)[2  +  i  -  noC"-,,--,)] 

+  aV(2  +  pJfV)  E-2 

llais  equation,  v/hich  takes  the  place  of  equation  l6  and  which  reduces 
to  equation  l6  when  h  (external)  =  0,  my  be  comnared  tern  by  term  v;ith 
equation  1^4.  In  equation  E-2,  quantities  have  the  meanings  ascribed 
to  them  on  page  5»  but  the  evaluations  are  at  the  outer  radius.  Mote 
that  all  sipxis  are  the  same  (as  in  equation  I'l)  excent  for  the  ^  terms. 
The  quantity  ^  is  the  external  aiTi;)ient  temperature  (for  convection). 

External  heat  transfer  can  take  place  both  by  convection  and  by 
radiation,  so  that  the  coeffcient  h  in  equation  E-1  is  the  sum  of 
two  coefficients 

h  =  h^  +  h^  E-3 

The  first,  h  ,  is  the  actual  convective  coefficient.  Vie  have  presumed 
c 

natural  convection  from  a  cylinder  into  still  air,  using  the  relation 

h  =  0.24t(e/d)^*^^  E-'4 

on  page  218  of  reference  8,  Here 

0  =  iTj^,  -  E-5 

and  d  is  the  diameter  of  the  cylinder,  l.e.,  the  outer  diameter  of  the 
pipe.  The  "radiative-convective"  coefficient  h^  is  an  equivalent  con¬ 
vective  coefficient  to  account  for  radiative  heat  transfer  and  is  de- 
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fined  by  the  equation 

Q  =  aeACCT^^^+O^-dJ^C)**]  =  h^ACT^-i^)  E-6 

in  which  aisthe  Stefan-Boltzmann  constant  (0. 171^1  xl0“®  Btu/hr  ft^  B**), 

G  is  the  emissivity  (we  take  O.R  for  bare  steel),  A  is  area  (’which  cancels), 
il)  is  the  teiTperature  of  the  surroundinr^^  with  which  radiative  heat  trans¬ 
fer  takes  place  (this  may  be  and  usually  is  different  from  the  convective 
ambient  temperature  i^),  and  C  =  ^6o  is  the  constant  required  to  convert 
Fahrenheit  teirperatures  to  Rankine  (absolute)  temperatures. 

■fhe  value  of  h  detemined  from  equation  E-3  is  in  units  of  feet  and 
hours  rather  than  inches  and  seconds  and  it  must  be  multiplied  by  1/(720)^ 
before  being  used  in  PIPETEIi.  There  is  provision  in  PIPETEfi  for  using 
the  relation 


h  =  f  h  +  f  h 
c  c  r  r 


E-7 


where  f^  and  f^  are  "augmenting  factors"  which  may  be  used  to  account 
for  the  fact  that  the  geometry  is  not  exactly  that  of  an  infinite  cylinder 
so  that  a  "face"  or  protuberances  such  as  nuts  and  bolts  mav  convect 
and/or  radiate  more  heat  than  would  otherwise  be  accounted  for. 

The  convective  ambient  temperature  ijj  and  the  radiative  ambient 
temperature  ^  could  be  made  to  be  functions  of  time,  .just  as  with  the 
internal  fluid  tenperature  (t>,  but  this  has  not  been  done  in  the  present 
inplementatlon  in  which  ip  and  ip  are  taken  to  be  constant. 

The  input  (and  the  input  instructions)  have  been  modified  slightlv 
to  accorrmodate  this  added  facility.  For  each  problem  an  additional  card 
is  inserted  after  the  first  card  described  previously.  If  this  card  is 
blank,  perfect  Insulation  is  assumed  at  the  outer  surface.  Other^'/ise, 
the  card  inputs  ip  ^  and  g. 
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Another  slif^t  modification  has  also  been  made  in  bIP^TI’’n.  Tn 
the  first  card  of  a  problem  decl<,  settinr  TIMIT  =  -500.  causes  the 
profTam  to  search  for  and  read  initial  temperature  data  as  described 
earlier,  while  setting  TINIT  to  any  value  greater  than  -^85.  caases 
the  program  to  assume  a  constant  initial  temnerature  eoaal  to  the 
value  of  TINIT.  The  new  capai)illty  is  that  setting  Til  ITT  =  -^100. 
causes  the  program  to  start  with  initial  temneratures  eoaal  to  the 
final  temperatures  for  the  preceding  problem.  Unless  the  values  of 
N  are  equal  for  the  two  problems,  some  nonsense  is  generateci.  In 
use,  of  course,  this  capability  is  intended  for  a  succession  of  tv/o 
or  more  problems  involving  the  same  pipe,  divided  into  the  same  num¬ 
ber  of  layers.  A  modified  set  of  COMTitllT  cards  at  the  beginning  of 
the  program  describes  the  modified  capabilities  and  the  details  of 
the  problem  deck. 
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